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Abstract 



A non-linear integral equation (NLIE) governing the finite size effects of excited states 
of even topological charge in the sine-Gordon (sG) / massive Thirring (mTh) field theory, 
deducible from a light-cone lattice formulation of the model, has been known for some 
, time. In this letter we conjecture an extension of this NLIE to states with odd topological 

charge, thus completing the spectrum of the theory. The scaling functions obtained as 
solutions to our conjectured NLIE are compared successfully with Truncated Conformal 
■ Space data and the construction is shown to be compatible with all other facts known 

, about the local Hilbert spaces of sG and mTh models. With the present results we have 

achieved a full control over the finite size behaviour of energy levels of sG /mTh theory. 



1 Introduction 

Finite size effects are widely recognized to play a major role in modern QFT (Quantum 
Field Theory) jjj]. The concept of scaling functions, encoding the dependence of energy levels 
on the finite volume (Casimir effect), has shown its practical convenience for a control of 
renormalization flows under change of scale of many QFTs. Most of the progress made so 
far in computing the scaling functions including non-perturbative effects have dealt with two- 
dimensional QFT, in particular with integrable QFTs, where exact methods a la Bethe Ansatz 
can be used to diagonalize Hamiltonians. A control of the behaviour of energy levels of a 
two-dimensional integrable QFT under change of the finite volume of a cylinder space-time 
geometry can lead not only to useful renormalization flow information, but also serves as a 
basis for future developments in the investigation of correlation functions @. 

A method of investigation that has been very fruitful, even for theories which are not 
integrable, is the Truncated Conformal Space (TCS) approach ||], which is an intrinsically 
non-perturbative method, but has problems of principal nature, coming from the fact that 
one does not have an analytic control of the spectrum, and of practical nature, because to 
reach a certain precision in the resulting energy levels one has sometimes to resort to very 
high truncation levels and introduce enormous matrices to diagonalize. 
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For integrable QFT, there also exist exact analytic methods to compute the finite size 
effects, like e.g. the Thermodynamic Bethe Ansatz (TBA) §. Another approach is the non- 
linear integral equation (NLIE) introduced as a continuum limit of Bethe equations emerging 
from a light-cone lattice regularisation, by Destri and De Vega |||, [|. (Similar equations had 
already appeared in a different condensed matter context in j?J). These methods provide an 
impressive numerical precision of the scaling functions, compared to the TCS one, and also 
give an analytic control. The principal limitation of these exact methods was basically that 
only the vacuum (and a few excited states degenerating with it at L — > oo) could be obtained. 
However, progress of the last few years have opened the investigation to excited states (see e.g. 
for the TBA approach refs. In the Destri-De Vega approach, the excited state equation 

was introduced in [ 10 1 and later generalized to include generic configurations of Bethe roots 
in 1 11]. The correct quantisation rule was established in [12, 13], where it was shown that the 
resulting scaling functions agree with numerical TCS data. 

The power of the NLIE method lies in the fact that it allows to describe all excited states 
in a single framework for all values of the sG coupling constant, and it can be used even in 
cases when the vacuum TBA is not known (in contrast with the TBA approach of [||, It 
was argued in |10] that a so-called a-twist following the ideas of Zamolodchikov |14] gives an 
access to scaling functions in $(1,3) perturbations of minimal models, the simplest examples 
of which include the scaling Yang-Lee and the Ising model in zero magnetic field. 

The NLIE still has a serious limitation, namely that only states lying in sectors with even 
topological charge can be accessed. This restriction is an apparently unavoidable feature of 
the light-cone lattice construction from which the NLIE is deduced. However, it is of primary 
importance to master also the scaling functions of the states with odd topological charge, as 
for example the single soliton state in sG, or the single fermion state in mTh. In this letter 
we give a conjecture of how the NLIE looks for such odd topological charge states, starting 
from the results for the UV conformal dimensions computed in [13]. The solutions of the 
NLIE proposed for the odd sector are then compared to the TCS approach for sG/mTh, i.e. 
perturbed c = 1 CFT that we have developed in |l^, |13"|| , Finally, the possibility to twist the 
NLIE itself leads us to a full control of all periodic/antiperiodic boundary conditions. This 
allows us to describe the differences between mTh and sG that should be expected in the odd 
charge sectors, in full accordance with the ideas of [15|. 



Throughout the paper we adhere to the same notations as in [12, The reader is invited 



to consult those papers (and references therein) as an introduction to the concepts and general 
setup. 



2 UV analysis in the NLIE framework 

The NLIE describing the excited state spectrum of sG/mTh theory for even value of the 
topological charge can be deduced from an inhomogeneous six-vertex model put on a light- 
cone lattice and has the following form in the continuum limit 

Z{0) = MLsmh& + g(0\d j ) + 2Sm I dxG(d - x - irj) log (l + {-if e iZ(x+i ^\ , (1) 

where L is the spatial volume and M. is the soliton mass. We also introduce the dimensionless 
volume parameter I = ML. The function g^l'&j) is the so-called source term, composed of 
the contributions from the holes, special objects and complex roots (close and wide) which we 
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call sources and denote their positions by the general symbol = {hk , Uk , Ck > Wk} and is 



N H N s M c M w 

k=l k=l k=l k=l 

where the function \ is the soliton-soliton phaseshift in sG theory, G is its derivative up to a 
factor of 27r, and the index II denotes the second determination necessary for wide roots (for 
details see fl3|]). The positions of the sources are related to the counting function Z({f) by the 
Bethe Ansatz quantisation conditions 
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where Ij are the Bethe quantum numbers. There is a so-called counting equation 

N H - 2N S = 2S + M C + 29(tt - 2 1 )M W , (2) 

relating the number of the different roots/holes to S, which is the spin of the XXZ chain 
appearing in the lattice regularisation and takes integer values for the Bethe states. The sG 
topological charge Q turns out to be Q = 25. 

Once the solution for the counting function Z{$) is known, the energy can be computed 
from the formula 

N H N S 

E = M V] cosh hj - 2M V] cosh yj 
j=i i=i 

M c M w 

(3) 



-M. ^2 cos ^ c j + -M. (cosl 1 w j )n 
i=i i=i 

/CO j r 
— 2Sm sinh(x + i77) log(l + (-l) s e iZ( - x+i ^) 
-oo 27r L 



with a similar expression for the momentum. 

We start our investigation by recalling some facts about the UV analysis of the spectrum 
of the NLIE, i.e. the consideration of the so called kink equation and its consequences. The 
UV limit of sG theory is c = 1 conformal field theory (CFT) with a compactification radius 
R which is related to the six-vertex anisotropy 7 by 



R- l = \h(i- x ) , (4) 



7T 



and it is often convenient to introduce the parameter^] 

7T 

P = 1 . 

7 

The primary fields are vertex operators Vr n m \ with conformal weights 

(n,m) 2 1 R 2 



4 p > 1 (i.e. < 7 < 7r/2) corresponds to the repulsive while p < 1 (i.e. 7v/2 < 7 < ir) to the attractive 
regime of sG theory. The nth breather threshold is at p — 1/n. 
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where m € Z is the winding number (which is conserved in the off-critical theory and is 
identified with the sG topological charge Q) and n is the field momentum. The relation ([|) 
means that we identify the perturbing operator with Vm ) + V(__ 10 )- 



We recall some of the results of ref. |13||. In the UV limit I — > the sources can be classified 
into three types: their position can remain finite ("central"), or they can move towards the two 
infinities as ±log| ("left /right movers"). We denote the number of right/left moving holes by 
N H and similarly we introduce the numbers N^, and M w . 
We introduce the partial spins S + , S~ by the definition: 



S ± = \[N±- 2iV± - M± - 2M±0(tt - 2 7 )] . 



The expression for the conformal weights is 

_ r± ^ jl 

2vr 16vr 2 (l -7/vr) ' 

where 



A± = ± (I± - 2I± - - 4) + g± + ^ , (5) 



cu ± = ±2(vr - 2 7 ) (5 - 25 ± ) - 2(vr - 7) (5 + 2fc±) , (6) 
fc± are integers determined by the condition — ir < uj± < it, 

S± = -45 ± (5-5 ± )f^^ + 27rg± , (7) 

1 — 7/7T 

and we introduced the following notation for the sums of quantum numbers of left/right 
moving sources of the different types 















T w = and if 






3=1 


3=1 





(q w is an integer or half-integer depending on the configuration of the wide roots, which is 
best calculated case by case). 

From (||) we obtain the following expressions: 



m = 2S 
2n = (5 + 2S + M sc ) mod 2 



(8) 



where M sc is the number of selfconjugate wide roots. Since S € Z we have access only to 
sectors with even value of the topological charge Q = m. In jDJ, the selection rule for sG 
states was shown to be 

(5 + M sc ) mod 2 = 0. 

Now let us consider putting S half-integer in (||) and remember that in order to have the 
sG with periodic boundary conditions we need n G Z. Thus we obtain the rule 



(5 + M sc ) mod 2 = 1. (9) 

From (§) we see that if S € Z + - then iV# must be odd (given that the number of close 
complex roots Mq is always even). Therefore to describe such states we put an odd number of 
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hole sources in the NLIE. This extension of the equation cannot be derived from the light-cone 
lattice approach, since the latter only allows for even number of holes. In the following we 
will take it as a conjecture and give arguments in support. 

Let us now present examples of UV dimensions computed from the NLIE for some simple 
states. For the state with one hole with quantum number 1 = 0, using 5 = 1 we obtain the 
conformal dimensions ^ 

A = 8^2 ' 

corresponding to the vertex operator V/Q iy This operator can be identified as the UV limit 
of the one-soliton field as it is well known [15]. 

Consider now the case of three holes, with quantum numbers I\ = 1, 1% = and I3 = —1. 
The conformal dimension we obtain corresponds to the vertex operator V(o,3) as expected. If 
we choose other values of quantum numbers different from this "minimal" choice, we obtain 
descendents of vertex operators Vr n ,3) , where n € Z depends on the number of holes moving 
to the left and to the right. 

As a last example, let us take three holes with quantum numbers I\ = 1, I2 = and 
I3 = —1, but we add a close pair with quantum number for both complex roots. We obtain 
the conformal weights 

a± = pp + 1 - 

which is the first nonchiral descendent of Vr 0> iy 

Using the above results, we formulate our conjecture as follows: the sG states with odd 
value of topological charge correspond to taking an odd number of holes and choosing the 
quantisation rule 5 in accordance with (||). In the next section we shall see that using this 
rule we get scaling functions which are in excellent agreement with numerical results obtained 
from the c = 1 Truncated Conformal Space (TCS) method. 



3 Numerical comparison between the NLIE and TCS 

In this section we present the numerical results which give further support to the conjecture 
on the identification of sG states with odd topological charge. In particular, we compare the 
NLIE prediction to numerical results obtained from the TCS method for the lowest energy 
levels. 

We start with the one-soliton state, which is conjectured to be a one-hole configuration 
with integer quantisation. Figure [I] shows the comparison graphically, while table [l] gives an 
idea about the numerical magnitude of the difference. Note that the deviations are extremely 
small for small values of I and they grow with the volume, exactly as expected for truncation 
errors. 

Next we take the two lowest energy levels corresponding to states of three solitons with 
the same charge (i.e. three holes with quantum numbers 0,±1 and 0, ±2, respectively). Here 
we only present the energy plot (figure ||) . The numerical deviations are exactly of the same 
magnitude as in the case of the one-soliton state. 

The value of the coupling constant at which the above data are presented lies in the 



attractive regime where the TCS method is free of UV divergences (12, |13| and therefore we 



can compare absolute energy levels, after subtracting the linear bulk energy term 

^ = Bl, 5 = -itan^, 
M ' 4 2' 



5 



6 - t 



Figure 1: Comparison of the results for the one-soliton state coming NLIE and TCS at p — =, The 
continuous line is the NLIE result, while the diamonds represent the TCS data. 
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NLIE 


TCS 


Relative deviation 


.5 


6.080571 


6.08062 


0.000008 


1 


3.126706 


3.12685 


0.00005 


1.5 


2.177411 


2.17791 


0.0002 


2 


1.727224 


1.72776 


0.0003 


2.5 


1.475004 


1.47593 


0.0006 


3 


1.320353 


1.32168 


0.001 


4 


1.153188 


1.15548 


0.002 


5 


1.075376 


1.07908 


0.003 



Table 1: Numerical comparison of the energy levels predicted by the NLIE to the TCS data for the 
case depicted in figure |]. The NLIE data are exact to the precision shown. 
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Figure 2: The energies as a function of the volume for first two three-soliton states coming from NLIE 
and TCS at p = |. The continuous lines are the NLIE result, while the diamonds represent the TCS 
data. 
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NLIE 


TCS 


Relative deviation 


.5 


25.27656 


25.3398 


0.0025 


1 


12.79207 


12.8854 


0.007 


1.5 


8.674948 


8.79114 


0.01 


2 


6.647771 


6.78091 


0.02 


2.5 


5.455906 


5.60033 


0.026 


3 


4.680766 


4.83188 


0.03 


4 


3.751729 


3.90497 


0.04 


5 


3.231486 


3.37598 


0.045 



Table 2: Numerical comparison of the energy levels predicted by the NLIE to the TCS data for the 
case depicted in figure ||. 



which is predicted by TBA QT6J| and also by the NLIE [§, from the TCS data. 

The state with three holes and one complex pair can be most easily investigated in the 
repulsive regime, since it changes its nature in the attractive regime where the complex roots 
behave in a more complicated way. The result is depicted in figure |3[ As in the repulsive 
regime the TCS is UV divergent, one can only compare relative energy levels, and as we are 
in the sector with unit topological charge, the natural way is to normalize the energy levels to 
the one-soliton state, which is the ground state of the Q = 1 sector. Although this cannot be 
directly seen from figure || it is clear from the numbers presented in table ^ that the truncation 
errors are much larger even after subtracting the UV divergences. The same phenomenon was 
observed in our earlier works [ 



4 Scaling functions in the massive Thirring model 

Here we formulate another conjecture concerning the finite size energy levels of the massive 
Thirring model. Of the two quantum numbers n, m specifying a vertex operator Vr n>m \, m can 
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Figure 3: The NLIE prediction for the state with three holes (with quantum numbers —1, and 0) 
and a complex pair, at p = 1.5, compared to TCS data. Both are normalized to the ground state of 
the Q = 1 sector, i.e. the one-soliton state, to get free of the UV divergent part of the TCS data. 



be identified with the sG topological charge and takes integer values, while n is not conserved 
in the massive theory. Generally, one can introduce twisted boundary conditions, for which 

where the parameter v lies between and 2ir. Such boundary conditions can be realised by 
adding a constant term a to the right hand side of the NLIE which is related to v as 

v = a — tt5 mod 2tt . 

Here we will need only the cases v = or tt, which are realised by taking n G Z/2. For 
convenience, we divide the Hilbert space into four sectors, given by 

ft/ = 0{^"(n,m) : n G Z, m G 21} , 
Hn = @{F M ■ n£Z, m G 2Z + 1} , 

Hni = ©{^(n,m) : n G Z + -, m G 2Z + 1} , 

f 

Hiv = 0{^"(n,m) : n G Z + -, m G 2Z} , 
where F( n ,m) denotes the Fock module corresponding to the vertex operator Vi n >m ) . According 



to |15(| , the sG Hilbert space (for periodic boundary conditions) can be identified with TLi®TLn 
while that of mTh is TLi®TLni- Note that these two spaces are related to each other by twisting 
the boundary conditions. The corresponding two algebras are two different maximal closed 
local operator algebras, which are essentially unique (see comments at the end of the section) 
at generic values of the compacitification radius R. Tiiv never appears as a local sector in 
sG/mTh theory. Using the results from the UV limit of the NLIE (|8|), we see that the massive 
Thirring model is selected by (putting a = 0) 

5 = M sc mod 2. (10) 
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1 


NLIE 


TCS 


Relative deviation 


3.0 


1.549935 


1.55081 


0.0006 


3.5 


1.415153 


1.41657 


0.001 


4.0 


1.321588 


1.32357 


0.0015 


4.5 


1.254597 


1.25722 


0.002 


5.0 


1.205454 


1.20881 


0.003 



Table 3: The comparison between the NLIE results for the energy of the one-fermion state and the 
TCS data for the range 3.0 < I < 5.0. The coupling constant is p — j and the hole quantum number 
was chosen I — \. 

As an example, let us consider the one-hole state with half-integer quantisation and quantum 
number I = — . Using the general formulae we compute the UV dimensions 

2 \2R 2 J ' 

which corresponds to the vertex operator Vn^i). If we choose instead I = — |, we obtain the 
UV limit V(_ 1/2,1). These are exactly the positive charge components of the (Dirac) fermion 
field in the c = 1 CFT and so we can identify the one-hole configuration with the one-fermion 
state of the massive Thirring model. One can also compute the dimensions of more general 
configurations and verify that they describe states in Hi © T~tni- 

The iterative solution of the NLIE presents serious difficulties when the arrangement of the 
source terms is not symmetric under i? — ► — The number of steps generally increases when I 
decreases. This can be prevented by making a re-antisymmetrization of the counting function 
computed at each step, which however can be done only for symmetric source configurations. 
For mTh states of odd charge, the source configurations are never symmetric. In the case 
of the one-fermion state, the iteration slows down too much to reach values of I less than 
3.0. The solution of this problem may necessitate the introduction of some fundamentally 
new numerical algorithm to solve the NLIE and prevents us from presenting a numerical 
comparison to TCS data of such mTh states for the time being. The results available for 
regions of I where the iteration converges fast enough (see table ||) show an agreement with 
TCS similar to the one presented for sG states in the previous section (cf. table [l]). 

Our selection rule (|l0|) takes into account that, as described in our previous paper, we 
redefined the source term for the self- conjugated roots (see section 3.2 of |jl3|). If one keeps 
the form of the source terms as they are derived from the lattice, one can summarize the 
selection rules for the even and the odd sector of mTh model in the following simple form: 

5 = mod 2 . 

With the same convention, the selection rule for the sine-Gordon model becomes 

5 = 25 mod 2 . 

Let us make a remark on the uniqueness of the maximal local operator algebra. The two 
operator algebras correspoding to the bosonic theory sG and the fermionic theory mTh 

A = {V(n,m) '■ n G Z, m G Z} 

Af = {VL m ) : n G Z, m G Z + n/2} 
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are the only ones at a generic value of the compactification radius R in the sense that any other 
maximal local operator algebras can be mapped into them by redefining R by multiplication 



with an integer |15]. However, the identification of the perturbation operator is not unique: it 
can be taken to be any of Vk,o + V-k,o with k integer. This would change the value of R in @ 
by a factor of k. If we map these theories back to perturbations by Vi t o + Vl^o by redefining 
R to the value in (||), we obtain the following other local operator algebras: 

A ( b k) = {V {n/Kkm) : meZ,neZ} 

-4/ fe) = {V(n/k,km) ■ m eZ, n eZ + m/2} 

This is the complete list of possible local theories as k runs over all positive integers. Note 
that sG/mTh corresponds to choosing k = 1. Another interesting theory is the one defined 

(2) 

by A b with the corresponding Hilbert space given by Hi © Ti.iv- At the Kosterlitz-Thouless 
point 7 = where the perturbation becomes marginal this describes the current-current 
perturbation of a level-1 SU{2) WZNW model, which is nothing else than the SU(2) Gross- 
Neveu theory. Note that the states with odd topological charge are absent from the local 
sector. This is due to the fact that the elementary excitation is a "kink" with fractional 
Lorentz spin 1/4. To describe the kink state, created by the operator V±xu±\ in the UV, one 
has to twist the odd sector NLIE by a = ±7r/2. This sort of nonlocality is a general property 



(k) 

of the models based on A b I which is why they were called "kink" theories in ||15|| . 

5 Conclusions 

In this letter we have proposed a generalization of the NLIE to the odd charge sectors of 
sG and mTh models. The conjecture has been tested against the data coming from a TCS 
diagonalization of the c = 1 perturbed CFT Hamiltonian that is expected to give a description 
of sG/mTh models. The agreement is excellent and gives strong evidence for our conjecture, 
which is also supported by the fact that the proposal is the most obvious and natural gener- 
alization of the even sector NLIE. Further evidence is given by the UV analysis reproducing 
correctly the conformal dimensions of the states to be found in such sectors, and in particular 
by the successful identification of the single sG soliton state as well as the two components of 
the mTh fermion. 

It is interesting to notice that the NLIE framework together with this proposal allows 
for a complete treatment of both sG and mTh Hilbert spaces. The two models, according 



to the interpretation of ref.[15|, correspond to two subspaces of the same space, identified in 
the UV limit by the action of two different maximal local subalgebra of conformal operators. 
One leads to a counting of states described by the modular invariant partition function, while 
the other corresponds to the T2-invariant one. mTh model corresponds to sG with twisted 
boundary conditions and vice- versa. All these features can now be controlled also away from 
UV criticality, thanks to the NLIE approach that we have in a sense "completed" with the 
present extension to odd charge sectors. 

Of course one could think to investigate different local UV operator algebras, identifying 

(k) 

other models, e.g. the algebras A b j. These are controllable on the level of the NLIE by 
choosing appropriately the twist parameter a implementing boundary conditions. Thus one 
could explore, e.g., the formulation of the deformed SU(2) Gross-Neveu model which is the 

(2) 

QFT having A b as local operator algebra away from the Kosterlitz-Thouless point. Other 



10 



interesting QFTs arise at special ("rational") values of the radii, like the $(1,3) perturbations of 



minimal models. Notice that in the language of ref. fli6 |, all these models are just "projections" 
onto different local sectors, of the same Hilbert space of sG model on the cylinder, if a-sectors 
are meant to be included in the Hilbert space from the beginning. 

It would be interesting, as a further development, to propose a modified light-cone lattice 
formulation of sG/mTh, leading to Bethe equations that allow the treatment of odd number 
of holes. A proposal in this respect has been put forward by Destri and Segalini |Q, in order 



to avoid the nonlocality present in the original light-cone lattice approach [18|. This nonlo- 
cality is manifested both in the lattice Hamiltonian and in the construction of the Thirring 
fermion field on the lattice. The fermionic light-cone lattice model has Bethe equations that 
allow for odd number of holes and look very similar to the ones derived from the "traditional" 
light-cone approach. However, since in this framework the fermion field is represented locally 
on the lattice, there appears the well-known doubling of fermionic species. The doublers are 
shown to decouple from the physical fields in the continuum limit of the Hamiltonian. In 
order to give a derivation of the NLIE starting from the fermionic formalism it is necessary 
to follow this decoupling in the Bethe equations themselves, which is out of the scope of this 
paper. We would just like to point out that the rule of half-integer Bethe quantisation in the 
mTh model is in complete accord with their fermionic lattice model. 
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